Delta Effects in Pion-Nucleon Scattering 
and the Strength of the Two-Pion-Exchange 
Three-Nucleon Interaction 



V.R. Pandharipande 1 , D.R. Phillips 2 , and U. van Kolck 3 



1 Department of Physics, University of Illinois, 
Urbana, IL 61801 

2 Department of Physics and Astronomy, Ohio University, 
Athens, OH 45701 

3 Department of Physics, University of Arizona, 
Tucson, AZ 85721 



Abstract 

We consider the relationship between P-wave ttN scattering and the strength of the 
P-wave two-pion-exchange three-nucleon interaction (TPE3NI). We explain why 
effective theories that do not contain the delta resonance as an explicit degree of 
freedom tend to overestimate the strength of the TPE3NI. The overestimation can 
be remedied by higher-order terms in these "delta-less" theories, but such terms 
are not yet included in state-of-the-art chiral EFT calculations of the nuclear force. 
This suggests that these calculations can only predict the strength of the TPE3NI 
to an accuracy of ±25%. 



PACS nos.: 21.30.Cb, 13.75.Cs, 12.39.Fe, 11.30.Rd 



1 Introduction 



A long-standing quest in hadronic physics is to relate the properties of free pions, observed 
in, for instance, pion-nucleon (ttN) scattering, to those of the pions which play such a 
significant role in the nuclear force. Recently, the Nijmegen group has provided a striking 
demonstration that one-pion exchange indeed provides the longest-range component of 
the two-nucleon potential. They extracted, with small error bars, the masses of the 
charged and neutral pions and the couplings of pions to the nucleon from fits to the pp 
and np scattering data [TJ. A subsequent Nijmegen analysis of NN data then confirmed 
that two-pion exchange [21 EJ EJ OH IB] gives a significant fraction of the intermediate-range 
attraction in the NN interaction [7j. (In some models other mechanisms, e.g. the very 
broad a meson [H], also contribute to this attraction.) In systems beyond A = 2 the 
three-nucleon interaction plays a subtle, but important, role. In this paper we focus on 
the Fujita-Miyazawa (FM) [§| term in the two-pion-exchange three-nucleon interaction 
(TPE3NI) . It appears — at least for light nuclei — that this is the largest piece of the three- 
nucleon force [TU] . 

Ideally nN scattering data should be used to directly construct the TPE3NI. However, 
the pions that generate nuclear forces are highly virtual. The relation between the scat- 
tering they experience from nucleons inside the nucleus and that observed in free space 
is non-trivial. To determine it, an extrapolation of the ttN amplitude from the "physical 
region" — where the pion energies are greater than m n — to the "virtual region" — where 
pion energies are much less than — is needed. 

The delta isobar is the most prominent feature of ttN dynamics. The delta peak in the 
7T + p elastic scattering cross-section is larger by an order of magnitude than any other jH] . 
Therefore, when constructing models of the ttN interaction that will be used for the 
extrapolation to the virtual region it is natural to include the delta as an explicit degree 
of freedom. This was the path followed many years ago, and the leading two-pion-exchange 
two- and three-nucleon potentials with an explicit delta were derived by Sugawara & von 
Hippel [2] and Fujita & Miyazawa (FM) jU], respectively. These two-pion-exchange NN 
and NNN potentials were recently re-derived as pieces of the more general expressions 
for two- and three-nucleon forces that are obtained when an effective field theory (EFT) 
with explicit delta degrees of freedom is applied to the problem of nuclear forces [3 HI] . 
Here we discuss how the FM potential arises in any theory with an explicit delta. Our 
expression for this potential is connected to ttN scattering data through the delta mass 
and the 7riVA coupling constant, both of which can be determined from the irN data. 

But the highly- virtual pions exchanged in the TPE3NI have energies much less than 
the delta-nucleon mass difference. This has encouraged the development of an approach to 
nuclear forces that is different from that of Sugawara & von Hippel and Fujita & Miyazawa. 
In this approach the delta degree of freedom — along with all other ttN resonances — is 
"integrated out". This yields an EFT in which pions and nucleons interact in the most 
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general way. In this EFT ttN interactions are point-like, and are organized as an expansion 
in the number of space and time derivatives (for a review, see Ref. |12j). The expansion 
parameter is essentially with u the pion energy and AM = — M pa 300 MeV 
~ Im^ the delta-nucleon mass difference. Applying this 'delta- less' EFT to 7rN scattering 
is challenging (see, e.g. Ref. [13]) since the expansion parameter is, at best, |, and the 
expansion breaks down completely at the delta peak. However the expansion should 
converge well if u <C AM, a condition which should have fair validity in nuclear-structure 
physics. The leading contributions to NN and NNN potentials in this EFT were found 
in Refs. [Hj and [TT] . respectively 1 . 

We have argued that nuclear-structure physics is within the domain of validity of both 
the theory with explicit deltas and the 'delta-less' EFT. We might expect then, that the 
two theories would give similar results for the strength of the TPE3NI. But this turns 
out not to be the case. Effective theories without an explicit delta predict a strength for 
the TPE3NI that is 1.5 to 2.5 times larger than that obtained by FM [T^]. Studies of the 
spectrum of light nuclei with the Green's function Monte Carlo method, including three- 
nucleon interactions, favor a strength of the TPE3NI closer to the FM value [TU1 |2*U] 2 . 

Here we identify the origin of this discrepancy. Parameters in the Lagrangian of the 
theory with pions and nucleons alone must be extracted from ttN scattering data. But 
the poor convergence of the derivative expansion in that theory tends to contaminate 
parameters extracted in this way. These parameters then appear in the TPE3NI and 
result in overestimation of its strength. Within the delta-less EFT this problem is only 
mitigated if many orders in the expansion are retained. 

This simple argument is presented as follows. In Section |21 we write down an EFT 
with nucleons, pions, and explicit deltas, and compute, to leading order, both the P-wave 
nN scattering amplitude and the TPE3NI. In Section El we use a theory without explicit 
deltas to compute the TPE3NI. By construction the ttN amplitudes in this theory and 
the theory of Section El agree at ttN threshold. We show that they differ by a factor of | 
in their prediction for the strength of the FM NNN potential. We then discuss how this 
overestimation would be remedied at higher orders in the delta-less EFT, and what the 
implications of this problem are for contemporary EFT computations of the TPE3NI. 

lr The delta contributions were of course implicit in previous dispersion-theoretical approaches |15lllfi| 
and models |17l I18| . although the correct chiral-symmetry properties are difficult to maintain when 
connecting the pion-nucleon amplitude to the potential without using field theory |19| . 

2 This conclusion is somewhat dependent on the regulator used in the three-nucleon force, but holds 
definitively if one requires that the cutoffs used in the NN and NNN system be the same. 
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2 A theory with explicit deltas 



Although many terms contribute to ttN scattering and the three-nucleon potential, here 
we focus on the delta contributions. We do not claim that this is an accurate or complete 
model for either ttN scattering or the TPE3NI, but it serves to illustrate the point we wish 
to make regarding the relationship between ttN data and the strength of the TPE3NI in 
delta-less EFTs. For discussions of this relationship in the context of hadronic models, 
see, e.g. Ref. [TH] . 

We consider P-wave ttN scattering in an effective theory with an explicit delta degree 
of freedom. We will be interested in small pion momenta, and so we need only the leading 
terms in the ttNN and nNA interaction Lagrangians. These are: 

CnNN = |fiVWiV-V$ (1) 

CnNA = |f (AtSTA + H.c.)- V$ (2) 

where N, and A are the pion, nucleon and delta fields, f n ~ 93 MeV is the pion 
decay constant, qa — 1.29 is the axial-vector constant that corresponds to the value of 
the (charged) nNN coupling constant reported in Ref. PQ, — 2.8 is the corresponding 
pion-nucleon-delta transition strength, and S and T are Rarita-Schwinger transition spin 
and isospin operators. Both S and T obey generalized Pauli identities of the form: 

S f A S B = -A B - -ia- A x B. (3) 

3 3 K ' 

Alternatively, one can work with the Hamiltonians 

H n NN = -^<T-V($(r)-T) , (4) 
H wNA = [S ■ V($(r) • T) + St • V($(r) • T+)] , (5) 



where, at this order, f nNN = and Lata = E rr A . 
2.1 ttN scattering at low energies 

At leading order in small momenta these Lagrangians yield four diagrams that contribute 
to P-wave irN scattering. They are shown in Fig. ^ Only two involve the delta. They 
give the nucleon-pole-subtracted amplitude that enters the TPE3NI. Graph A.l is the 
direct — or s-channel — graph, and graph A. 2 is the crossed — or w-channel — graph. 

We evaluate these graphs in the center-of-mass (COM) frame in which the pion energy 
is u, and denote the momentum and isospin of the initial (final) pion by qi and ti (q2 and 
t 2 ). Since we limit ourselves to pion momenta of the order of the pion mass the nucleon 
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Figure 1: Four irN scattering diagrams. Dashed lines represent pions, solid line nucleons, 
and thick solid lines delta isobars. 

kinetic energies are smaller than w by a factor of order rn^/M, and can be neglected in 
this leading-order calculation. For the same reason we neglect the kinetic energy of the 
delta. 

The delta contribution to the ttN amplitude is then given by 
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The Xj an d Xj are spin-isospin quantum numbers of the nucleon before and after scatter- 
ing. 

Using Eq. (jSJ), we can rewrite this amplitude as 
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2.2 The three-nucleon scattering amplitude 

We now turn our attention to the tree-level delta contribution in the TPE3NI. To this 
end we consider the amplitude for nucleon % emitting or absorbing a pion of momentum 
±qi and isospin ti and nucleon k emitting or absorbing a pion of momentum ±q 2 and 
isospin t 2 . In "direct" diagrams the pion "1" converts nucleon j to a A and "2" reconverts 
it to nucleon. In the "crossed" diagrams "2" converts and "1" reconverts. There are 12 
"direct" and 12 "crossed" diagrams in time-ordered perturbation theory. The 12 direct 
diagrams are shown in Fig. 

The contribution of the direct diagrams to the three-nucleon scattering amplitude is 
given by: 

f 2 / - \ r 12 
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Figure 2: Twelve "direct" NNN diagrams. Notation as in Fig. ^ 

Here Xij,k an d Xi,j,k denote the initial and final spin-isospin states of nucleons i,j and 
k, and U a is the product of the three energy denominators in diagram a of Fig. |21 The 
values of n a can be read off the diagrams, and they are listed in Table Once again we 
have neglected nucleon and A kinetic energies in computing these denominators, which 
is valid in our leading-order calculation. 
From Table ^we can easily verify that: 
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Substituting this in Eq. (jSJ) gives: 
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Table 1: The values of (— 1/U a ) for direct diagrams. 



The contribution of the crossed diagrams involves analogous energy denominators, and 
can be calculated similarly. The sum of direct and crossed diagrams, 

P ( -1 \ 

A ™ = -^Ffafckfc • q 2 T fc • t 2 \ Xk )(Xi\<7i ■ qiTf ■ tilxi) ( ^2AM J 

x^(xJ|S] • 02 Sj ■ qi T] • t 2 T, ■ t! + S] • q 1 S i ■ q 2 Tj ■ tiT; ■ t 2 | X ,),(ll) 
gives the Fujita-Miyazawa potential V^ k ' FM jH], 

t r2ir,FM fnNN f ^ \ j. j. 

x { J 1ftm) ( qi " q2ti ■ t2 - \°' • qi x q2Tj ■ ti x t2 ) - (12) 

This result agrees with many previous re-derivations of the FM potential, e.g. Ref. ^T]. It 
is exact at tree level in the static limit if the only terms in the irNN and 7riVA Lagrangians 
are those in Eqs. (JH) and (J2J). 



3 Relation to theories without explicit deltas 

We now attempt to find a more direct connection between ttN scattering data and 
V^' — one that does not invoke the delta as an explicit degree of freedom. Such at- 
tempts have been reviewed in Ref. whose notation we follow below. 

A key aspect of this connection is that nN scattering involves pions with u ~ m^, 
while in V^ k ' FM we have u ~ m^/M. (The typical nucleon momentum in the nucleus is 
of order the pion mass, and the pion energy is then smaller by a factor m n /M.) Since we 
have already been neglecting terms suppressed by m^/M we take q\ = q 2 = 0. Given this 
kinematics, the three-nucleon potential can be written: 

t>2tt _ fnNN a > " Qi " 02 



ai3 



(13) 
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where Ui = yqf + m^. comes from the pion propagators and 

- J** = 5^ [a + b qi • q 2 + c(g 2 + g 2 2 )] - d^e ^ ■ qi x q 2 ) (14) 

is the Born-subtracted nN subamplitude. The first term is due to S-wave scattering, the 
second gives the anticommutator part of the TPE3NI, the third is zero, and the fourth 
gives the commutator part. The first term is very small in the context of V//1 [TO], and it 
is zero in the present model. 

The crucial point, then, is the determination of the coefficients b and d. In a theory 
without explicit delta fields, they are fitted to irN data near threshold. If we lived in 
a world where there were no contributions to ttN scattering other than from the s- and 
ti-channel delta and nucleon poles, comparing Eq. (fTlj) and Eq. (JJJ) shows that a fit to 
threshold nN data would result in 



b = 4d =-£*±±( 2AM ) 
ml 9\(AM) 2 -mlJ 



(15) 



The TPE3NI corresponding to this amplitude is given by: 



V$ = ^J- 2 <T t -C ll <T k . q2 T l -t 1 T k . t 2 Of . (16) 



1^2 



The factor \ju)\uj\ comes from the pion propagators, and the factors besides OJ describe 
the coupling of the pions to the nucleons i and k. The nN interaction is described by: 

0] N = b fqi • q 2 ti • t 2 - i (Tj • qi x q 2 Tj ■ t x x t 2 J , (17) 

with b given by Eq. (|15|) . Of course, this is just the usual FM form, but with specific 
choices for the coefficients b and d. 



3.1 The problem 

Comparing the V?l in Eq. (jTfij) with the the "exact" result for our model (V^ of 
Eq. (|12|)) we find that they are the same apart from the crucial fact that the strength of 
the interaction in the "delta- less" theory has the factor r^^z^p > instead of the of 
the "exact" result. Since AM ~ 2m n , these factors are ~ and ~ respectively. 
One way to understand this result is to realize that the direct term for the 7rN scattering 
amplitude in Eq. (JJJ) and Fig. [T] is evaluated at the energy of a real pion, and so has the 
energy denominator AM — m n for low-momentum pions. This denominator is half of 
the average denominator, AM, of the diagrams in Fig. [21 that contribute to the TPE3NI. 
The crossed pion term mitigates this discrepancy, but not enough to cure the problem. 
Ultimately, the that is extracted "directly" from ttN scattering data is too strong by 
a factor of 4/3. 
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The difference between and V*l' FM is of order (j!^) 2 - It will vanish in the limit 
AM ^> m^, which includes the chiral limit m n — ► 0. However, in the context of the 
nuclear many-body problem m n is not small. The range of OPEP is comparable to 
the mean inter-nucleon spacing in nuclei, and the energies required to excite nucleons to 
isobar states such as the delta are not much larger than m n . 

Of course, in the real world there are contributions to the ttN amplitude other than the 
two graphs we have considered here. Also b and d will probably be determined from data 
that are not exactly at threshold. While we cannot say a priori in which direction these 
effects go, fitting ttN data at higher energies will presumably only make the extrapolation 
problem worse. 

Parts of this problem have been understood for a long time, but, as discussed in the 
introduction, the prevailing folklore has been that an EFT without explicit deltas could 
still work well in nuclei, because the relevant energies in nuclear-structure physics are 
much smaller than AM. However, the poor convergence of the EFT without explicit 
deltas for irN scattering affects the TPE3NI because b and d are not calculated from first 
principles; instead they are fitted to threshold ttN data. This necessitates an extrapolation 
from pion energies u ~ m n to the energies of the highly-virtual pions in the TPE3NI, 

2 

which are of order This extrapolation takes place over an energy range that is sizable 
compared to the radius of convergence of the "delta-less" theory — AM. 

Here we have explicitly considered the implications of such an extrapolation for the 
three-nucleon potential, but other few-nucleon potentials (including the two-nucleon force) 
will be afflicted by the same problem. All use 7rN parameters that are potentially con- 
taminated in a similar way. Such contamination will occur in all EFTs for low-energy 
hadronic physics which contain only pion and nucleon degrees of freedom. 

3.2 The solution 

In a theory with explicit deltas this extrapolation is under much better control, since the 
pion-energy dependence of the nN amplitude is better reproduced. In contrast, at leading 
order in the "delta-less" theory the coefficients of the two operators in Of are energy 
independent, and so the value extracted for them at threshold, where ui = m T , is used in 
the TPE3NI, where oo ~ 0. 

But at higher orders in this EFT additional corrections to the ttN amplitude, and in 
particular to the two operators in Of , enter. To see what form this higher-order energy 
dependence would take, we expand the result (JJJ) in powers of (^7) 2 • The first correction 
to the leading-order results for b and d (fTKJ) occurs 

at °[(s5?) 1- The form of °j N 

is now: 

Of = (b + L> 2 ) K • qjj ti • t 2 - ~ (Tj ■q 1 xq 2 Tj ■ t x x t 2 j . (18) 
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In the EFT, terms such as but 2 and duo 2 appear in the Lagrangian as pion-nucleon in- 
teractions with time derivatives. We must fit tcN data over a range of pion energies to 
determine both b and b. If, once again, we imagine living in a world where the true answer 
was given by Eq. (J7J), then fitting the form (fTSj) to reproduce (|7J) in the region around 
uo = m n yields: 

Af 2 NA ( 2AM \ / ml \ 

!~ / a STTvT —o ) i ( 19 ) 

(20) 



9 ml \(AM) 2 - ml J \ (AM) 2 - ml 
4 2AM 
9 m 2 ((AM) 2 -m 2 ) 2 ' 



Note that at ttN threshold this gives exactly the same result for 0^ N as in Eq. (|17|). 
However, extrapolating to = now yields a TPE3NI that has an additional factor of 

2 

(1 ~~ (AMp-m 2 ) * n ^ s strength. If we set AM = 2m 7r , this gives an overall factor of 
instead of the — found in the "exact" calculation with explicit deltas. This means that 
in the theory without explicit deltas the "exact" factor — is being built up as: 

J. = (x * (21) 

m w \ 6 9 / 3m„- 

a series that converges moderately quickly. 

To summarize: in the theory without explicit deltas it is important to realize that the 
factor ^p- obtained by fitting nN "data" with the leading-order form (|T7|) is not the final 
answer. This result will change when higher-order terms are incorporated in the theory 
and used to improve the extrapolation from u ~ m^ to to ~ 0. We can estimate the 
size of such terms based on our knowledge that the convergence will be governed by the 
parameter JnFr, and that — due to crossing — only even terms in this expansion can appear 
in b and d. The leading-order result should therefore be quoted as: 



i _ _ 4 fwNA 



9 ml 3m„ 



(22) 



More conservative error bars are certainly acceptable, but the ~ 25% we have chosen is 
the minimum permissible theoretical error that can be assigned to b when it is extracted 
in the theory without explicit deltas. Such an error bar turns out to be consistent with 
the "exact" answer for b in the simple model considered here. 



4 Conclusion 

We have shown that theories without an explicit delta tend to overestimate the delta 
contribution to the TPE3NI. This is because there is an error in the leading-order com- 
putation of the three-nucleon potential in the "delta-less" theory. The error is ~ 25%, 
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and it is necessary to include terms suppressed by (^j) i n the EFT to reduce it. The 
inclusion of other higher-order effects, such as nucleon recoil and dispersive effects for 
intermediate-state deltas, may make the extrapolation error smaller than we found, but 
it seems unlikely that it will completely remove the difficulty 

Unfortunately this problem is present in the state-of-the-art N 3 LO chiral EFT compu- 
tation of NN and NNN potentials [21]. The terms that ameliorate the overestimation 
appear in C^/ N , and so will not enter the chiral EFT nuclear force until N 4 LO. Comput- 
ing the two- and three-nucleon potentials to this (or higher) order will take considerable 
effort. It may well be that an EFT with explicit deltas is simply a more efficient tool 
than one without. In fact, the first studies in nuclear EFT [HI E] included diagrams 
with intermediate deltas in their calculation of the nuclear force. The drawback of such 
a treatment is that in order to fix parameters one must analyze data around the delta 
resonance, which necessitates a resummation of the delta self-energy. Only recently has 
a power counting been devised that allows a systematic EFT treatment of effects in this 
kinematic region [22J. 

The delta-less EFT has also found difficulties with certain nN parameters that are large 
because the effects of the integrated-out delta are encoded there. In Ref. j2H] Epelbaum et 
al. argued that there is a cancellation of delta-excitation and 7rp-exchange contributions 
in nuclear forces. This motivated their use of NN and NNN potentials containing ttN- 
interaction parameters smaller than those extracted from chiral analyses of 7rN scattering 
data. We stress that the reduction in strength of the NNN force we have discussed here 
is not based on such an argument. It is independent of details of nuclear dynamics at the 
distance scale l/m p . 

So, until the theory with explicit delta degrees of freedom is further developed, or 
delta-less theories can be extended to higher order, the irN parameters used in the NNN 
potential should be viewed as only loosely constrained by ttN data. Furthermore, EFT 
extractions of ttN parameters from NN data (see, e.g. Ref. [Zj) and from ttN data (see, 
e.g. Ref. JT3J) can be expected to give results that differ by amounts of order • 
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